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We calculate numerically universal finite-size-scaling functions of the magnetization for the three-dimensional 
0(4) and 0(2) spin models. The approach of these functions to the infinite-volume scaling functions is studied 
in detail on the critical and pseudocritical lines. For this purpose we determine the pseudocritical line in two 
different ways. We find that the asymptotic form of the finite-size-scaling functions is already reached at small 
values of the scaling variable. A comparison with QCD lattice data for two flavours of staggered fermions shows 
a similar finite-size behaviour which is compatible with that of the spin models. 



1. INTRODUCTION 

Our motivation for studying 3d O(N) spin 
models is their relation to quantum chromody- 
namics (QCD). The QCD chiral phase transition 
for two light-quark flavors is supposed to be of 
second order in the continuum limit and to be in 
the same universality class as the 3d 0(4) model 
jlj. In the staggered formulation of QCD on the 
lattice a part of chiral symmetry is remaining and 
that is 0(2). The comparison of QCD lattice data 
to 0(4) scaling functions for V — > 00 and expo- 
nents yields good agreement for Wilson fermions 
0, but not for staggered fermions ||. This may 
be caused by using too small lattices in QCD, 
so we will compare in this paper staggered QCD 
data to the finite-size-scaling functions of the 3d 
0(2) and 0(4) spin models. 

The 0(iV)-invariant nonlinear er-model on a d— 
dimensional hypercubic lattice is defined by 

/3H = -J J2 s * ,s i - H E Si ' W 

<i,j> i 

Si is an N-component unit vector at site i with 
a longitudinal (parallel to the magnetic field H) 
and a transverse component 

S, = gfft + Sj; . (2) 
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The order parameter of the system, the magneti- 
zation M , is given by 

M = <^£ 5 f> = <s">> (3) 

i 

and the longitudinal susceptibility by 

XL = m = v{{sf) ~ M2) - (4) 

2. FINITE-SIZE-SCALING FUNCTIONS 

The general form of the finite-size-scaling func- 
tion for the magnetization is given by 

M = L- & l v ^(tL x l v , hl}l v %L- w ) , (5) 

where t and h are the normalized reduced tem- 
perature t = (T — T c )/Tq and magnetic field 
h = H/Hq. On lines of fixed z = th" 1 ^ 5 and 
after expanding $ in L -UJ we have 

M = L- fs ' u Q z {hL 1 /^) + ... , (6) 

where Q z is a universal function. Examples of 
lines of fixed z are the critical line (z = 0) and 
the pseudocritical line (z — z p ), the line of max- 
imum positions of the susceptibility xl m the 
(t, /i)-plane for V —> 00. 

For large lattices Q z will converge to the 
asymptotic form Q Zt00 

Q z ^Q z ,oo = faiz^hL 1 ^) 1 / 5 . (7) 
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Figure 1. Lines of fixed z = 0.5, 1, 1.5, 2 (dashes), 
the pseudocritcal line (solid) at z p = 1.33 ± 0.05 
and from measured peak positions (squares). 

Qz,<x> can be calculated from the scaling function 
Jg for V — » oo, which we have calculated in Ref. 
{§ for 0(4) and in Ref. § for 0(2). 

3. NUMERICAL RESULTS 

Our simulations were done on three-dimensio- 
nal lattices with periodic boundary conditions 
and linear extensions L — 8 — 120 using the clus- 
ter algorithm. In order to eliminate corrections 
to scaling we used increasingly larger values of L. 
For the critical coupling and the exponents we 
took the same values as in Refs. § and @. 

First we had to determine the value z p of the 
pseudocritical line. The scaling function of \l 
for V — > oo can be determined from the scaling 
function of the magnetization fa 

Evidently, the maximum of \l at fixed h and 
varying t is at the maximum point z p of the 
function in the brackets of Eq. (||). We obtain 
z p = 1.56(10) for 0(2) and z p = 1.33(5) for 0(4). 

As a check we have also determined the peak 
positions on lattices with L = 24 — 96 which ex- 
trapolate for L — > oo to z p = 1.65(10) for 0(2) 
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Figure 2. (a) Finite-size scaling of ML^^ for 
0(4) on the pseudocritical line. The solid line 
shows Q Zp ,oo, the symbols denote different lattice 
sizes L. (b) is a double-log plot of (a). 



and z p = 1.35(10) for 0(4). In Fig. [I] we com- 
pare both methods for the 0(4) model. We see 
that both methods agree very well for small h, 
but for larger h the extrapolated peak positions 
yield slightly higher pseudocritical temperatures. 

In Fig. ||a we show the finite-size-scaling plot 
for 0(4) on the pseudocritical line (at z p = 1.33). 
There are strong corrections to scaling and the 
universal function Q Zp is approached from above. 
If one looks at the logarithmic plot (Fig. ||b), one 
finds that Q Zp converges to Q Zp ,oo from below 
and coincides with the asymptotic form at about 
HL 1 ^" « 20. 

We have also investigated the finite-size-scaling 
functions of 0(2) on the pseudocritical line and 
of both models on the critical line. They have a 
very similar behaviour to the one shown above, 
but on the critical line the corrections to scaling 
are negative, and in the case of 0(4) they are 
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Figure 3. versus m a aL 1 ^-, with 0(4) 

exponents, from QCD lattice date with two stag- 
gered fermions. The open (filled) symbols are 
data from Ref. § (@, multiplicated by 2). The 
lines are drawn to guide the eye. 

much smaller than on the pseudocritical line. The 
asymptotic form is always approached at small 
values of HL 1 ^" from below. 

4. COMPARISON TO N f = 2 QCD 

In Fig. m we show a finite-size-scaling plot of 
staggered QCD data on the pseudocritical line 
from the JLQCD collaboration || and the Biele- 
feld group 0. The exponents are from the 0(4)- 
model. Here the quark mass m q plays the role 
of the symmetry-breaking field h, and the chi- 
ral condensate {iJjiJj) is the order parameter for 
uiq — » 0. (iptp} has been taken at the peak posi- 
tions of the chiral suceptibility Xm a t fixed m q for 
each lattice size. The error from the peak posi- 
tion determination is not included in the plot. As 
in Fig. ||a the universal finite-size-scaling function 
is approached from above. At small m q the data 
are however higher than expected. 

In the logarithmic plot, Fig. [|, we compare the 
QCD data to a line c+ {1/5) h\(m q aL 1 / Uc ), which 
represents the asymptotic finite-size-scaling func- 
tion. Because of the unknown normalization of 
QCD c was chosen freely. There is good agree- 
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Figure 4. Logarithmic plot of Fig. |[ The solid 
line is the asymptotic finite-size-scaling function 
of 0(4). 

mcnt of the QCD data with the expected be- 
haviour, especially when one takes into account 
that the lattice sizes are still small and correc- 
tions are probably present. 

The corresponding result for 0(2) is very sim- 
ilar to that of 0(4) and, because of the spread of 
the QCD data, one cannot really distinguish the 
two cases. 

An extended version of this paper can be found 
in Ref. g. 
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